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Abstract 

We consider a variational problem with a polyconvex integrand and 
nonstandard boundary conditions that can be treated as minimization of 
the strain energy during the suturing process in the plastic surgery. Ex¬ 
istence of minimizers is proved as well as necessary optimality conditions 
are discussed. 
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1 Introduction 


Given an open bounded connected domain C C with a sufficiently regular 
(locally Lipschitz) boundary, dC, let us consider the integral 

I {u) := I W {S/u{x)) dx (1) 

Jn 

to be minimized on a class of Sobolev functions u : ^ ^ with a kind of 
boundary conditions to be described later. All over the paper we assume the 
integrand W : — i- RU {+oo} to be polyconvex. This means that the 

representation 

w^(G=5(T(G), 

holds for some convex function g : R'^G’^) —>■ Ru {-boo}, 


dAN 

T{d,N) :='^x{s), k{s) 



dim 

{s\)^{d-s)\{N-s)V 
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where 


T iO ■■= (AdjiC, Adj^C, AdjgC,..., Adj.^^O , ? e 

and Adjj,^ is the vector of all minors of the matrix ^ of order k = 1,2,... ,dAN, 
respectively. In particular, Adj^^ = ^ and Adj^;^ = det ^ whenever d = N. 

It is known that, under strong coercivity assumptions on W to assure weak 
convergence of the minors of gradients for the minimizing sequence, the func¬ 
tional I attains its minimum on m (•) -|- p > 1. We refer to the 

fundamental work by J. Ball [3] motivated by problems coming from nonlinear 
elasticity and to [DiiniiiH] for further improvements. 

The lower semicontinuity for general polyconvex integrands with respect 
to the weak convergence in IT^’^(r2; R.'^), ft C has been the subject of 
many investigations. Namely, Marcellini showed in |19| that this property holds 
whenever p < N. Later, this result was improved by Dacorogna and Marcellini 
in [S] who proved the lower semicontinuity for p > TV — 1 while Maly in [TH] 
exhibited a counterexample for p < iV — I. The limit case p = N — 1 was 
adressed in [D [H H HU. Very recently (see m and m) the limit case was 
studied for polyconvex integrands depending on x and/or on u. 

Besides the Dirichlet boundary condition u = u for the displacement one 
considered, for instance, boundary condition on traction, which somehow de¬ 
pends on the normal derivatives of u(-) (generalizing the Neumann boundary 
data). Observe that these conditions (displacement, traction or a mixed one) 
can be applied either to the whole boundary dfl, or to some of its subsets of pos¬ 
itive Hausdorff measure leaving the rest free. Moreover, some restrictions on the 
jacobian may be relevant and practically justified. For example, the constraint 
det Vu (x) > 0 means that the minimum is searched among the deformations 
preserving orientation, while det Vu (x) = 1 refers to the case of incompressible 
elastic body. 

One of the possible applications of the above variational problem is regarded 
to plastic surgery, namely, in the woman breast reduction, where we deal with 
a sort of very elastic and soft tissue. Some recent publications (see, e.g., [H 
m [211 US m) were devoted to mathematical setting of the related problems 
and to their numerical simulations. Medical examinations allow to consider the 
involved tissue as a neo-Hookean compressible material (see [23] )■ We have a 
more precise model when the strain energy is defined by the integral o with 
the density W : —>■ R, 

W (0 := P (tr (C • C"’) - 3 - 2 In (det o) 

-bA (det ^ - I)^ -b /3 tr (^Adj f • Adj , (2) 

where ”tr” means the trace of a matrix, Adj^ := Adj 2 ^, and the symbol ”T” 
stands for the matrix transposition. One of the steps of the (breast reduction) 
surgery is the suturing, which mathematically can be seen as an identification 
of points of some surface piece r+ C dfi with points of another one T” C dft. 
Denoting the respective correspondence between the points of r+ and T" by cr. 
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we are led to a new type of constraint 

u{x) = {uo a)(x), X € F'*", (3) 

called the knitting boundary condition. Let us note that the one-to-one mapping 
a is not a priori given and should be chosen to guarantee the minimum value 
to the functional ©■ In other words, a minimizer of ([T]) (if any) should be a 
pair (m, cr) where u G W^’P(n;IR^), p > 1, and a : F'*' —^ F is sufficiently 
regular. We set the natural hypothesis that a and its inverse a~^ are Lipschitz 
transformations (with the same Lipschitz constant L > 0). Practically this 
means that the sutured tissue can not be extended nor compressed too much. 

Motivated by the problem coming from the plastic surgery we will consider 
just the case p = 2 and d = N = 3, although the results remain true in the case 
p > 2 and arbitrary d = N > 2 as well. 

The paper is organized as follows. In the next section we give the exact 
setting of the variational problem together with the main hypotheses on the 
integrand W. For simplicity of references we put here also some important 
facts regarded with the Sobolev functions. In Section 3 we justify first the 
well-posedness of the problem by showing that the composed function from the 
knitting condition ([3]) belongs to the respective Lebesgue class. Afterwards, we 
prove existence of a minimizer as an accumulation point of an arbitrary mini¬ 
mizing sequence (the so called direct method, see [5]). The paper is concluded 
with a necessary optimality condition for the given problem (see Section 4) al¬ 
lowing to construct effective numerical algorithms, which can be successively 
applied in the medical practice. 


2 Main hypotheses and auxiliary results 


In what follows we fix a nonempty open bounded and connected set II C 
whose boundary dkt is assumed to be locally Lipschitz (see, e.g., [T71 p. 354]). 
By the symbol (dx) we denote the Lebesgue measure in the space K™, 
m = 2,3, while T-L^ means the two-dimensional Hausdorff measure (see, e.g., 

USD- 

Let us divide the surface dkl into several parts F^, i = 1,2, 3,4, in such a 
way that (Fi n Fj) = 0 for * ^ j. Moreover, we set F 4 := F+ U F“ where 
F=*= C F with V? (F=*=) > 0 and V? (F+ n F“) =0 are also given. 

Suppose that W : —>■ R is a poly convex function satisfying the growth 

assumption: 

bF(0 >co + ci|a"+C2|AdjfD+C3(detO", ^ G ( 4 ) 

where Cq G R and > 0, i = 1, 2,3, are some given constants. Here and in what 
follows by I • I we denote the norm of both a vector in R" and a 3 x 3-matrix. 


3 


Taking into account that tr j = |^|^ for each matrix ^ S we see 

that the integrand (I2|) satisfies the above properties. Indeed, it is convex as a 
function of T (^) being represented as a sum of three terms, which are convex 
w.r.t. det^ and Adj^, respectively. Furthermore, 

W{0 = -3fi + X + fi lei" + /3 |Adj ei' + /(det C), C £ 


where the function 

f (t) :=^t^ — 2Xt — 2^\nt, t > 0, 

is lower bounded by some (negative) constant. 

Since on various pieces of the surface dfl the boundary conditions are struc¬ 
turally different (some part of dfl can be left even free), to set the problem we 
use the notion of the trace operator, which associates to each u £ (n;]R^) 

a function Tru defined on the boundary, dXl, which can be interpreted as the 
” boundary values’’’’ of u. We refer to m pp. 465-474], where the existence 
and uniqueness of the trace operator were proved for scalar Sobolev functions 
u £ 'W^’P{n),p > 1. For vector-valued functions u : ^ instead, we can 

argue componentwise. So, applying m Theorem 15.23], we define the trace as 
the linear and bounded operator Tr : (n;R^) ^ L^(i9r2;]R^) satisfying the 

following properties: 


1 . 

2 . 


TrM(a;) = u{x), x £ dfl, whenever u £ fl C 


for each u £ W^’^(n;R^) and any test function (p £ 0^(12; R^) the equal¬ 
ities 


[ Uj-^dx = - f ipj^^dx+ f (pjTr{uj)i^id'H‘^, 

Jn oxi Jq oxi Jqq 

hold for each i,j = 1, 2,3 where ly := (lyi, ^ 3 )^ means the unit outward 
normal to dXl. 


In addition to the properties above, observe that the trace operator gives a 
compact embedding into the space (9fl;R^) that will be crucial to obtain 
the main result in Section 3. Namely, the following proposition takes place. 

Proposition 1 Let 12 C R^ &e an open bounded set with locally Lipschitz bound¬ 
ary. Then for each {un} C (S2;R^) converging to u weakly in (12;R^) 
the sequence of traces {Tru„} converges to Trn strongly in (512; R^). 

The proof is based essentially on the following lemma giving a nice estimate 
for the surface integral of the trace operator. 

Lemma 1 Let 12 C R^ &e as in Proposition [H Then there exists a constant 
C > 0 such that 

f iTrul'^d-H'^ <C (- [ \ufdx + e [ iVnl^dx) (5) 

JdQ Jn Jn / 

for any e > 0 and any u £ (12; R^) . 
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Proof. Given x € dil, due to the Lipschitz hypothesis there exists a neighbor¬ 
hood Ux of X such that Ux H dil can be represented as the graph of a Lipschitz 
function w.r.t. some (local) coordinates. Without loss of generality, we may 
assume that 


Uxr]n = {(y', j/s) : fx W) <Vs< fx iy') + 5x, y' G Gx} 


and 

Uxndn = {iy'Jx{y')):y' gGx}, 

where Sx > 0, Gx C is an open set in the space of the first two coordinates 
and fx ■ Gx —> M is a Lipschitz function. By compactness there exists a finite 
number of points a;* G dfl , i = 1,... ,q, such that 


= U {Uxi n dn). 


i=l 


Set St := Sxi, Ui := Uxi, Gi := Gxi and fi := fxi, i = l,-.-,g. Denote by 
L > 0 the biggest Lipschitz constant of the functions fi. 

Let us choose e < min{dj : i = 1,..., y} and consider first the function u G 
(D) n (D; K.^) . Given i G {1, ..., y} , by the Newton-Leibniz formula, 
for each x' G Gi and X 3 G K with fi {x') <xz<fi {x') -I- e we have 

du 

u{x',fi{x'))=u{x',X3)- / -— {x',s) ds. 

Jfi{x') 0x3 


In turn, by Cauchy-Schwartz inequality. 


\u{x,fi{x'))\ < 2 \u{x',X3)\ + 


/ rfi{x')+e 

du , , 

V) 

V/iG') 




< 


( 2 
2 |M(a;',a;3)| +e 

V "'bG') 




Hence, 


< 


|u (x', Mx'))f \/l + \Vft (x')l^ 

2VT+L^(\u{x',X3)\^ + e 

V JfiG') 


du , 



Integrating both parts of the previous inequality in the cylinder 

fit := {{x , X3) : / {x) < X3 < f {x') + e, x' & Gi} CUiC] ft, 


gives 


£ 


'c/inan 


\u (x)l^ (x) 


< 2a/1-KL2 



|m (x)!^ dx -I- £^ 





< 


2\/1 + L2 



\u {x)\^ dx + £^ 
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Summing in i = 1,..., g we have 


< 



(6) 


The inequality ([5]) for each u € follows now from ([6]) by the prop¬ 

erties of the trace operator and by the density of smooth functions. ■ 

Proof of Proposition 1. Given a sequence {un} C converging 

weakly to u in there exists M > 0 with 


I {x)\^ dx < M 


for all n G N. Then, by the Rellich-Kondrachov theorem (see [171 Theorem 11.21, 
p. 326]), Un ^ u strongly in (fl; . Applying Lemma[T]to Tr (u„ — u) = 
Tr Un — Tr u, we have 



ufdx + £ IVWr; 

JQ 


u\‘^dx + AeM). 


\7u\^dx) 


Hence 

limsnp j |TrM„ — Tm]^ d'h? < AeCM. 
n^oo J dQ 

Letting e —>• 0+ concludes the proof. ■ 

We will use also the so called generalized Poincare inequality (see [SJ Theorem 

6.1-8, p. 281]). 

Proposition 2 Given an open bounded domain H C with locally Lipschitz 
boundary and a measurable subset T C dfl with (T) > 0 there exists a con¬ 
stant C > 0 such that 


f \u (cc) ^ dx <C 

/ 1 Vu (a;)]^ dx -f 

j Tr u (x) dT-L^ {x) 

2' 

J 

Q. 

J 

n 

J 

T 



for eachueW^’^ 


Let us formulate now the boundary conditions in terms of the trace operator. 
Consider first a surface 5 C defined by some continuous function h : R^ —?> R, 

S := {u G : h (u) = 0} . 
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Then, given T > 1 we denote by El (r+; F ) the set of all functions a : r+ F 
satisfying the inequalities 

^\x-y\<\a{x) - a{y)\<L\x-y\ (8) 

for all x,y € F+, and introduce the set Wl C (0;]R^) x El (F+; F“) of all 
pairs (u, a) satisfying the (boundary) conditions: 

(Cl) Tru (a;) = a; for "H^-a.e. a; G Fi; 

(C 2 ) /i(TrM(a:)) = 0 for "H^-a.e. x G F 2 ; 

(C 3 ) Tru (a;) = Tru (ct (a;)) for H^-a.e. a; G F+. 

Thus, we can write the knitting variational problem in the form 

min|y W {Vu{x))dx : {u,a) . (9) 



r2 m(F2), h{x) = 0 


Figure 1: General scheme of plastic surgery. 

Let us emphasize the difference between the boundary conditions (Ci) —(C 3 ) 
above (see Fig. 1). The first condition (Ci) means that the surface Fi is a part 
of the elastic body (breast) that remains fixed. The condition (C 2 ), instead, 
can be interpreted as the knitting of a part of the incised breast (surface F 2 ) to 
the fixed surface S of the woman’s chest, while (C 3 ) means the knitting of the 
cut breast surface F 4 := F+ U F“, of F+ into F“. Finally, the piece F 3 of the 
boundary dQ remains free and admits an arbitrary configuration depending on 
the knitting process. 


3 Existence of minimizers 


Before proving the main existence theorem let us justify that for each a G 
El (F+;F~) the boundary condition (C 3 ) makes sense. 
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Lemma 2 Let C be an open bounded connected domain with locally Lip- 
schitz boundary and C Oil be T-L^-measurable sets with TL^ (F^) > 0 and 
7^^(F“''nF“) = 0. Then for each u £ and a £ Si;,(F“'';F“), 

L > 1, the composed function Tru o a : F'*' —>■ is measurable w.r.t. the 

measure TL^ on F+. 

Proof. Given u £ (fi; R^) by the density argument there exists a sequence 

of continuous functions Un : 17 —>■ R^ converging to u in the space (17; R^). 

Consequently (see the properties 1 and 2 of traces), Vn = Trun —)• Tru, as 
n —>■ 00 , in 1? (917; R^). Then, up to a subsequence, we have 

Vniy) ^Tru{y) Vy £ F“ \ (10) 

where C F“ is some set with null Hausdorff measure. So, it remains to 
prove that 

(u-i (E^)) = 0 (11) 

because in such case we deduce from cni) that the sequence of (continuous) 
functions {u„ {a (a;))} converges to Tr u {a (x)) for all x £ F+ up to the negligible 
set of points Eq := a~^ i^o)- 

On the other hand, uni) follows easily from the left inequality in (jS]) and 
from the definition of Hausdorff measure (see [131 P- 171]): 

CXD 

(E) := ^ linr inf ^ (diam Ai)^ 

i=l 

where the infimum is taken over all coverings of E with the diameters 

diam^i < e. In fact, given t] > 0 we find e > 0 and a family with 

Ufci O Ef , diam^i < sfL and 

OO 

^(diamAi)^ < 

Since due to dH) obviously diam (ct“^ (^i)) < Ldiam < e, i = 1, 2,...; 
the family (A^)}^^ is a covering of Eq and 

CXD OO 

(diam (tr”^ (A^)))^ < (diam A^)^ < rj, 

i=l i=l 

we conclude that i^o) = 0; Ih® 77^-measurability of x i->- Tru(cr(x)) 
on F+ follows. ■ 

Let us give now an a priori estimate for the ’’weighted” integral 

J |TrM(cr (x))|^d77^ (x), 

r+ 

implying, in particular, that the composed function Tr m o ct belongs to the class 
L2 (9f7;R3). 
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Lemma 3 Let f2 C and C dfl be such as in Lemma[H Then given L >1 
there exists a constant > 0 such that the inequality 



|Trw (cr {x))\^ dTL^ (x) < £,l j 


\Tru{y)\^dn^ {y) 


holds whenever u G and cr G Ei (F''";r ). 


( 12 ) 


Proof. To prove the estimate m we employ the local lipschitzianity of the 
surfaces r=*=. Namely, given y G F" we choose an (open) ball B {y,ey), Sy > 0, 
such that F“ni? (y, Sy) can be represented as the graph of a Lipschitz continuous 
function with respect to some (local) coordinates. Without loss of generality 
we can suppose that this function (say /“) is defined on an open set D~ from 
the cartesian product of the first two components x' := {xi,X 2 ) and admits as 
values the component 0 : 3 , i.e., 


F- n i? (y, Sy) = { {z', f- (/)) :z'eD-}. (13) 

By the compactness of F“ one can find a finite number of points y^,..., y'J G F“ 
with 

F- = F-n (14) 

Set Sj := Syj , DJ := D- and /“ (z') := f~ (z'), z' G DJ , j = I,... ,q. 

Similarly, for any x G F+ there exist > 0, an open domain C and 

a Lipschitz function /+ : £)+ —>■ K such that 

r+riB {x, 6,) = {{z', /+ {z')) -.z'eDt}. 

We do not loose generality assuming that the value of /+ is the last component 
of the vector 2 ; G F+ (as in O). Again due to the compactness argument there 
exists a finite number of points x^,... ,x'^ G F"*" such that 


F+ = F+n U B{x\6,) 

i=l 


where 

:= imm^ 1^) , i = l,...,r. (15) 

Set also := HA and (z') := f^i {z'), z' G Df, i = 1,..., r, and denote by 
Lr the maximal Lipschitz constant of the functions f~ : DJ —> K., j = 1,..., y, 
and /+ Df -)• R, f = 1,..., r. 

We claim that given i = 1,..., r and cr G the image a (F+ n B (a:*, ^i)) 
is contained in F“ ni? (y^, Sj) for some j = 1,..., y. Indeed, let us choose j such 
that cr (a:*) € B (y^,Sj/2) (see (fT4l) l. Taking an arbitrary z G F+nB (a:*, (5i) we, 
in particular, have \z — x*| < |j; (see (fT51) l and by ([5]) |cr (z) — a (a;*)| < ej/2. 
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However, assuming that a[z) f^B{%j\ej) we have |cr( 2 :) —y^\ > Sj and 

hence 

\a{z)-a{x^) \ > \a{z)-y^ \ - |cr(x*) - 2 /^ | > Ej - y = y, 

which is a contradiction. In what follows we associate to each a G Sl and to 
each z = 1 ,..., r an index j = j (a, i) G { 1 ,..., g} such that 

a (r+ n B {x\Si)) CT~ r\B {y\ej) . 

The latter inclusion allows us to define correctly the (injective) mapping 
: Df —> D~, such that 

'I'cr * j ( ct . i ) 

cr(a:) =a{x'J+ (x')) = (i}\{x'))) ,x' G D+. (16) 

^From (|5]) it follows immediately that is Lipschitz: 

\'<P\{x')-'>P"a[z')\ < \a {x', f+{x')) - a {z', {z'))\ 

< L {\x'- zf + {f+{x') - f+{z')Yy (17) 

< L-iJ 1 + Lp \x' — z '\, x',z'GDf. 

So, by Rademacher’s Theorem, tplr is differentiable on with /Im¬ 

measurable gradient, and the inequality 

(x')| <M := L^l + Ll 


holds for a.e. x' G Df. Notice that the inverse mapping (z/;),) ^ is well defined 
on Gl, := tpl, (Z3)*') C by the formula similar to (fTHll . namely, 

(y')) = ((V'a) 

In the same way as Czl we deduce, from ([8|), that (V'^) is Lipschitz on the 
(open) set G), and, the estimate 


\y') 


< M 


(18) 


holds for C?-&.e. y' gG]^. 

Integrating the function |Tru (cr (a;))|^ on the surface piece r+ n H 
we pass first to the double integral 


J \Tiu {a {x))f dH^ {x) 

r+nB(x\Si) 

= Jj |Tru (cr {x'Jt (x')))|^ \/l+ |V/+ {x')\^dx' (19) 

D+ 

< v^l + Lp JJ jTru (cr (x , {x')))\'^ dx'. 

D+ 
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Due to the representation (ITBl) we can make the change of variables y' = (^0? 

x' S Df, in the integral (IT^ . and returning then to the surface integral on a 
piece of F”, we have 



Tr u (cr (a;',/+ {a:')))|^ dx' 


< 


< 


JJ \TTu{y'Jj {y'))\^ dety (iIjI) ^ (y') dy' 

Gi 

6M^ JJ\Tru {y\ /" (y')) |' ^I + fj {y')\^dy' 

Gi 

6M^J \TTu{y)fdH^ (y). 


r- 


( 20 ) 


Here we used the estimate (fT^ and the obvious inequality |detH| < 6 ||A||^ {A 
is an arbitrary 3 x 3-matrix). Since the sets r+ fl H (a;*, , i = 1 ,..., r, cover 

the surface r+, taking into account (IT^ and (l 20 l) we conclude that 



iTru (cr {x))\^ dTL^ (x) 


< J \Tru{a{x))\^ dn^ (x) 

’'^^r+nB(x\Si) 

< £j \TTu{y)fdn^y) 


r- 


where £ := 6rM^^l -|- Lp > 0 depends just on the Lipschitz constant L > 1 
and on the properties of the domain H (namely, of its boundary). ■ 


Proving the existence theorem we pay the main attention to the validity of 
the boundary condition (C 3 ) where Lemma |3] is crucial. 

Theorem 1 Let W : ^ RU {-boo} be a polyconvex function satisfying the 

growth assumption 0- Then problem m admits a minimizer whenever there 
exists at least one pair to := (tt, cr) G Wl with 

I (uj) := J W (Vu(x)) dx < -boo. 
n 

Proof. Let us consider a minimizing sequence {{unjCTn)} C Wl of the func¬ 
tional o, e.g., such as 


J W (Vun (x)) dx < inf 
n 


W (Vu(a:)) dx : {u,a) G Wl 


1 

-I— < -boo. 
n 


( 21 ) 
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Taking into account the estimate (HI) we deduce from m that the sequences 
{Vu„}, {Adj Vu„} and {det Vm„} are bounded in (fi; and in (fi; R), 

respectively. Applying Proposition [5] and the boundary condition (Ci) we find 
a constant C > 0 such that the inequality 


J \un (a;)|^ dx < C 


\VUn ix)f dx + 


X dH^ {x) 


holds for each n > 1. So, the sequence {un} is bounded in (ri;R^) and 
by the Banach-Alaoglu theorem, up to a subsequence, converges weakly to some 
function u € (r2;R^). Without loss of generality, we can also assume that 

{Adj Vu„} and {det Vu„} converge weakly to some functions ^ S (r2;R^^^) 
and 77 € L^(f];R), respectively. Now, by Theorem 8.20 [H pp. 395-396] due 
to the uniqueness of the limit we deduce that ^ {x) — Adj Vm (x) and rj {x) = 
detVu(a;) for almost all x G fl. Thus we have the weak convergence of the 
sequence {T (Vu„)} to the vector-function T (Vu) in the space 

On the other hand, recalling that {cr„} C Sl (r''';r“) (see ([5])) by Ascoli’s 
theorem up to a subsequence, not relabeled, {cr„} converges uniformly to ct G 

El (r+;r-). 

Since the integrand W is polyconvex, it can be represented as W (^) = 
g (T (^)), ^ G R^^^, with some convex function g : R, and, therefore, 


JW (Vm (x)) dx = 
n 

< 

Thus, it remains just to prove that oj := (u, ct) G Wl (he., that the Sobolev 
function u satisfies the boundary conditions (Ci) —(C 3 ) above with the trans¬ 
formation a). The validity of (Ci) and (C 2 ) follows immediately from Propo¬ 
sition [T] In fact, the weak convergence of {un} in (n;R^) implies the 
strong convergence of traces {Tru„} in (911; R^). So, up to a subsequence, 
TrM„(x) —>■ Tr'i((x) for H^-a.e. x G dfl. In particular, TrM(x) = x and 
h{Ti'u{x)) — 0 almost everywhere on Pi and on r 2 , respectively (w.r.t. the 
Hausdorff measure). 

In order to verify the condition (C 3 ) we observe first that 

Tvunix) = TrM„(cr„(x)), n = I, 2 ,..., ( 22 ) 

for x G r+ and consider the surface integral 

J := [ jTr-it (x) — Tru (ct (x))|^ (x). 

Jr+ 

By the Minkowski’s inequality we have 

jl/2 < ^jny/2 ^ ^jny/2 ^ ^jn^l/2 ^ 23 ) 


/ 5 (T (Vu (x))) dx < liminf / g (T (Vm„ (x))) dx 
J J 

n o 

inf I (^)) S / • 
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where 


J"” := [ \Tru{x)-TTUn{(Jn{x))f d'H'^ (x); 

Jr+ 

J 2 ■■= [ |TrM„ (cr„ (cc)) - Tru(cr„ (x))|^ (cc); 

Jr+ 

■= [ \TTu{an{x))—Tru^d {x))f dT-L^ (x). 

Jr+ 

Taking into account the equalities (1^^ by using Proposition [T] we immedi¬ 
ately obtain that J7i" —>■ 0 as n —>■ 00 . 

Due to the linearity of the trace operator, applying Lemma [3] and again 
Proposition [T] we arrive at 

1/2" ^ ■Cl y |Tr (un — u) {x)f dH^ (a:) —>• 0 as n —>■ 00. 

Let us approximate now u G W^’^(r2;R^) by a sequence of continuous 
functions : D fc = 1 , 2 ,... (with respect to the norm of (D; K^)). 

Then (see Proposition m Vk = TiVk ^ Tru as /c —>■ 00 in (i9D;R^). In 
particular, given e > 0 there exists an index fc* > 1 such that 

[ |TrM(y) - Ufc. (y)l^d'H^ (y) < e. 
jr+ 

By using Lemma [3] similarly as was done to estimate the integral J 2 have 


[ \TTu{anix)) - Vk- {(7n{x))f d'H'^ {x) 

Jr+ 

< \Tvu{y)-Vk-{y)\^d%^{y)<ilLe, n = l,2,..., (24) 


and similarly 



|Tru (tr (x)) 


Vk- (cr (x))!^ dT-L^ (x) < £,lS. 


(25) 


On the other hand, by the uniform continuity of Vk- and the uniform convergence 
(T„ —>■ CT as n —>■ 00 , we find a number n* > 1 such that 


\vk- {un {x)) - Vk- (cr (a;))| < e 
for all n > n* and all x € P'*", and, consequently, 

/ jvk-(cr„ (x)) - Vk-(d-(x))l^ d?{^ (x) < ?{^ (r~'') £^,n > n*. (26) 

Jr+ 

Joining together the inequalities (I24L (I33|) and (Ei) we obtain that 

< (Cl£)^/" + (Cle)^^" + (r+) n*. 

Since £ > 0 is arbitrary and the constant £l does not depend on n = 1,2,..., 
we conclude that all the three integrals in the right-hand side of (l23l) tend to 
zero as n —)■ 00 . Thus JT = 0, or, in other words, Tr'iZ (x) — Tru (ct (a;)) = 0 for 
Td^-a.e. X G P’*', and the theorem is proved. ■ 
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4 Necessary conditions of optimality 


In this section, under some additional hypotheses, we deduce necessary condi¬ 
tions of optimality for problem (|^. 

To simplify, assume that the function W is twice continuously differen¬ 
tiable and h is continuously differentiable. Moreover, suppose that the surfaces 
Ti,r 2 ,r+,T”,r 4 are sufficiently smooth. 

Given T C dfi, with 'H^(r) > 0, in what follows we denote by C^(r+;R^) 
the family of restrictions to T of all functions it : 17 —)• R^, whose gradients 
are continuous up to the boundary. Let us supply C^(r7;R^) with the natural 
sup-norm. 

We consider the problem dH) defined in the space C^(f7;R^) x C^(r+;R^). 

Theorem 2 Let (UjCT) € C^(r7;R^) x C^(r+;R^) be a minimizer of problem 
m- Assume that Vh{u{x)) ^ Q, x € T 2 and detVit((T(a;)) 7 ^ 0, x € 17. Then 


the following conditions are satisfied: 

Div(VW)(ViZ(x)) = 0, xef7; (27) 

VW(ViZ(x))t^(x) = 0, X e Lg; (28) 

VlT(Vu(x))ti(x) X Vh{u{x)) = 0, X e Lz; (29) 

VlT(Vu(x))ti(x) = 0, xer±. (30) 


Proof. Let us write the constraints in the minimization problem ([^ as F{u, a) = 
0 where the map F : C\n-,R^) x Ci(r+;R3) ^ Ci(ri;R3) x Ci(r 2 ;M) x 
C^(r+;R^) is given by 

F{u, a) := (it(x) — x, h(u(x)), u(x) — it(cr(x))). 


DF{u,a){u,d-){x) = 


Under our assumptions the map F and the functional I are both Frechet differ¬ 
entiable. In particular, for the (Frechet) derivative of F at the point (it, a) we 
have 

u{x) 

(Vh(it(x)), u{x)) 
u{x) — u(d'(x)) — Vu(a(x))d'(x) 

Here, and in what follows, (•,•) denotes the inner product in R^. Taking into 
account that V/i(it(x)) 7 ^ 0 on F 2 and that the Jacobian matrix Vit(d’(x)) is 
not degenerated, we have that the linear operator DF{u,a) is onto the space 
C^(ri;R^) X C^(r2;R) x C^(r+;R^). By the Lagrange multipliers rule (see, 
e.g., m) there exist linear continuous functionals Ai, A 2 , A’*’ on C^(Fi; 1 
C^(r 2 ;R) and C^(r+;R^), respectively, such that 
» 3 


E 


dW(yu{x)) dui{x) 




dxj 


dx+Xi{u)+\ 2 {{Vh{u),u))+X^{u—u{cr)—Vu{a)a) = 0. 


Applying the Divergence theorem we get 


(VIF(Vm(x))ii(x),m(x)) dLL^^x) 


— / (Div(VIF(Vi((x))), ■u(x)) dx-I- 
Jn Jan 

+Ai(i() - 1 - A 2 ((Vh(u), u)) + A^(u — u{a) — Vu((t)ct)) = 0. 


(31) 
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Here v{x) is the unit outer normal to the boundary. Varying u in (IHTl) such 
that u{x) = 0 on 9f2, we obtain (l27l) . Taking then u G with u{x) = 0 

on dVL \ Fa we arrive at (1281) . Furthermore, choosing appropriate functions u in 
m we obtain 


lr2 


(yWiyu{x))v{x),u{x)) dV?{x) + X 2 {{'Vh{u),u)) = 0 


(32) 


whenever u{x) = 0, x £ Oil \ r 2 , and 

f {\/W{yu{x))i'{x),u{x)) (^^{x) + — u{d)—Vu{d)(T) = 0 (33) 

dr+ur- 

whenever u{x) =0, x £ dQ \ (r+ U F”). 

Denote by F^ the part of r 2 where the vectors a{x) := VW{Vu{x))i^{x) and 
b{x) := V/i(u)(x) are co-linear. Taking an arbitrary c £ C^(r 2 ;M) such that 
c{x) = 0 and Vc(x) = 0 on F 2 , let us define 


u{x) 


a(x){a(x),b(x))-b(x)\a(x)\'^ 
(a(a:),6(ai))^ —|a(a;) |^|6(ai)|^ ^ ^ ’ 

0, 


a^er2\ro, 
X £ r§. 


Obviously, ii £ C^(F 2 ;K^), (u(x), a(x)} = 0 and ('u(x),b(x)) = c(x) for x £ r 2 . 
From (15^ we get A 2 (c) = A 2 ((&,u)) = 0. Hence, varying u £ C^(r 2 ;R^) in (15^ 
in a suitable way (in particular, setting u(x) = 0 on F®) we get a(x) = 0 in 
r 2 \ F 2 . Thus the equality (1^ follows. 

Finally, taking u = 0, from (l33l) we get A'*' = 0, and, as a consequence. 



(VW(Vu(x))iy(x),u(x)} d7{^(x) + [ (VkF(Vu(x))i/(x),-S(x)) (i'H^(x) = 0, 

Jr+ 


which implies (15(11) . ■ 
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